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MODELLING WOLBACHIA INFECTION IN A SEX-STRUCTURED 
MOSQUITO POPULATION CARRYING WEST NILE VIRUS 

JOZSEF Z. FARKAS, STEPHEN A. GOURLEY, RONGSONG LIU, AND ABDUL-AZIZ YAKUBU 


Abstract. Wolbachia is possibly the most studied reproductive parasite of arthro¬ 
pod species. It appears to be a promising candidate for biocontrol of some mosquito 
borne diseases. We begin by developing a sex-structured model for a Wolbachia in¬ 
fected mosquito population. Our model incorporates the key effects of Wolbachia 
infection including cytoplasmic incompatibility and male killing. We also allow the 
possibility of reduced reproductive output, incomplete maternal transmission, and 
different mortality rates for uninfected/infected male/female individuals. We study 
the existence and local stability of equilibria, including the biologically relevant and 
interesting boundary equilibria. For some biologically relevant parameter regimes 
there may be multiple coexistence steady states including, very importantly, a co¬ 
existence steady state in which Wolbachia infected individuals dominate. We also 
extend the model to incorporate West Nile virus (WNv) dynamics, using an SEI 
modelling approach. Recent evidence suggests that a particular strain of Wolbachia 
infection significantly reduces WNv replication in Aedes aegypti. We model this via 
increased time spent in the WNv-exposed compartment for Wolbachia infected female 
mosquitoes. A basic reproduction number Rq is computed for the WNv infection. 
Our results suggest that, if the mosquito population consists mainly of Wolbachia in¬ 
fected individuals, WNv eradication is likely if WNv replication in Wolbachia infected 
individuals is sufficiently reduced. 


1. Introduction 

Wolbachia is a maternally transmitted intracellular symbiont, and it is the most com¬ 
mon reproductive parasite infecting a significant proportion of insect species, see e.g. 
[IIIIM]- Wolbachia typically inhibits testes and ovaries of its host, and it is also present 
in its host’s eggs. It interferes with its host’s reproductive mechanism in a remarkable 
fashion. This allows Wolbachia to successfully establish itself in a number of arthropod 
species. Well-known effects of Wolbachia infections include cytoplasmic incompatibility 
(Cl for short) and feminization of genetic males also known as male killing (MK for 
short), see e.g. 0 [121 [31] |32]■ Another important well-known effect of Wolbachia infec¬ 
tions is the inducement of parthenogenesis, see e.g. ISED]. All of these contribute to the 
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fact that the mathematical modelling of Wolbachia infection dynamics is both interesting 
and challenging. 

In recent decades a substantial number of mathematical modelling approaches have 
been applied to model different types of Wolbachia infections in a variety of arthropod 
species. Perhaps most frequently researchers have been focusing on the development of 
mathematical models for Wolbachia infections in mosquito species. Many of the earlier 
models took the form of discrete time matrix models, written for population frequencies, 
see e.g. [331131], and the references therein. Using frequency-type models a number of 
researchers investigated the possibility of coexistence of multiple Wolbachia strains, each 
of which exhibits different types of the reproductive mechanisms mentioned earlier, see 
e.g. miiniiiniiMj- Among others, Wolbachia strains have been investigated as a poten¬ 
tial biological control tool to eradicate mosquito borne diseases. Originally the focus has 
been on Wolbachia strains that induce life-shortening of their hosts. This is because for 
many vector borne diseases only older mosquitoes are of interest from the point of view 
of disease transmission. Therefore the use of (discrete) age-structured population models 
has become increasingly prevalent, see e.g. |28j and the references therein. Fairly re¬ 
cently, in [33] the results of laboratory experiments were reported envisaging a successful 
introduction of a life-shortening Wolbachia strain in the mosquito species Aedes aegypti. 
In [33] three key factors, namely, strong Cl, low fitness cost and high maternal trans¬ 
mission rate, were identified as drivers of a successful introduction of the new Wolbachia 
strain into an Aedes population. To this end researchers have developed and analysed 
continuous age-structured population models for Wolbachia infection dynamics, which 
take the form of partial differential equations, see cni; which can often be recast as delay 
equations, see e.g. [Hdl]. 

In recent years there have been substantial modelling efforts to theoretically investi¬ 
gate the potential of biological control tools for limiting the impact of mosquito borne 
diseases. It is now widely recognised that biological control represents a viable alter¬ 
native to established methods such as the use of insecticides and bed nets. Among 
others, the sterile insect technique has been investigated in the recent papers [3 [31 [33]. 
More recently, it was reported that particular strains of Wolbachia (completely or almost 
completely) block dengue virus replication inside the mosquito hosts, see for example 
mmm- To this end Hughes and Britton [19] developed a mathematical model for 
Wolbachia infection as a potential control tool for dengue fever. Their work suggests that 
Wolbachia may be effective as such a control measure in areas where the basic repro¬ 
duction number Rq is not too large. These recent results underpin the possibility that 
Wolbachia may be a promising candidate for biocontrol of mosquito borne diseases, in 
general. Besides dengue. West Nile virus (WNv) is another well-known mosquito borne 
disease of current interest. WNv infection cycles between mosquitoes (especially Culex 
species) and a number of species, particularly birds. Some infected birds develop high 
levels of virus in their bloodstream and mosquitoes can become infected by biting these 
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infectious birds. After about a week, infected mosquitoes can transmit the virus to sus¬ 
ceptible birds. Mosquitoes infected with West Nile virus also bite and infect people, 
horses, and other mammals. However, humans, horses, and other mammals are ‘dead 
end’ hosts. This virus was first isolated in the West Nile region of Uganda, and since 
then has spread rapidly, for example in North America during the past 12 years. Since 
there is no vaccine available, the emphasis has been mainly on controlling the vector 
mosquito species. Some recent experiments, see m, have confirmed that replication of 
the virus in orally fed mosquitoes was largely inhibited in the wMelPop strain of Wol- 
bachia. Interestingly, in a recent paper, Dodson et al demonstrated in laboratory 
experiments that the wAlbB Wolbachia strain in fact enhances WNv infection rates in 
the mosquito species Culex transalis. However, in [5] the Wolbachia was not a stable 
maternally inherited infection, but rather they infected transiently somatic mosquito tis¬ 
sues, and hence the wAlbB infection did not induce significant immune response in the 
mosquitoes. This is probably key to their findings. Here we will focus on modelling a 
maternally inherited Wolbachia infection in a population model, which hypothesizes a 
large number of successive generations. Nevertheless, the findings in [B] underpin the 
importance of Wolbachia research in general and highlight the importance of contrasting 
the findings of new theoretical, laboratory and field investigations. 

In this work we introduce sex-structured models for Wolbachia infection dynamics 
in a mosquito population. This will allow us to incorporate and study the well-known 
effects of Cl and MK of particular Wolbachia infections, simultaneously. First we will 
treat a model which only involves the mosquito population itself. Then we will use this 
model as a basis for a much more complex scenario incorporating WNv dynamics in a 
Wolbachia infected mosquito population. The full WNv model will naturally include the 
bird population, too. 


2. Model for a Wolbachia infected mosquito population without WNv 

2.1. Model derivation. We start by introducing a model for a Wolbachia infection in 
a sex-structured mosquito population, incorporating sex-structure using a well estab¬ 
lished approach originally due to Kendall [2T]. More recent papers of Hadeler [12] and 
Hadeler et al [T3] derive and discuss sex-structured pair formation models in depth. We 
only model (explicitly) the adult population of mosquitoes. Our model allows us to take 
into account the well-known effects of cytoplasmic incompatibility (Cl), incomplete ma¬ 
ternal transmission, fertility cost of the Wolbachia infection to reproductive output, and 
male killing (MK), at the same time. We note that it was shown in [9] that a stable coex¬ 
istence of MK and Cl inducing Wolbachia strains is possible, in principle. Introduction 
of male killing Wolbachia strains in vector populations may have a significant effect on 
the disease dynamics, as typically only female mosquitoes are transmitting the disease. 
Also note that according to [d5| , those Wolbachia strains which cause greater disruption, 
as in the case of dengue transmission, confer greater fitness costs to the mosquitoes. This 
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may well be the case for West Nile virus, hence we account for the reduced reproductive 
output in our model. 

We deduce our starting model from basic principles. In particular, first we deduce 
mating rules arising at the individual level. Starting with the adult population of size iV, 
we construct a random mating graph. This is a bipartite graph, not necessarily complete, 
in which each vertex has degree at most one. The vertices represent male and female 
individuals and edges represent realized matings. Let us denote by M,Myj,F,Fjjj the 
numbers of un/infected males and females, respectively. For every adult mating pair, 
offspring is created according to the following rules. Below, m, f and denote 

uninfected/infected male/female individual, respectively. The parameter /? models the 
reduced reproductive output of Wolbachia infected females, r measures maternal trans¬ 
mission in the sense that it is the probability that a Wolbachia infected mother passes on 
the infection to its offspring, q measures Cl in the sense that when a Wolbachia infected 
male mates with an uninfected female, q is the probability that there is no viable off¬ 
spring. Finally, 7 measures MK in the sense that it is the probability that a Wolbachia 
infected male larva dies during its development. A complete list of parameter values will 
be given later on. With this notation the mating rules are described below. 


(1) m X f : create one pair of the same type (to, /). 

(2) m X fw ■ with probability /3, create no offspring. This reflects the fecundity 
reduction due to the Wolbachia infection. In the complementary case, with 
probability (1 — /3 )t(I — 7 ), create a new pair (to^j, /m), at the same time with 
probability (1 — create (0,/„,), i.e. a female only brood. This accounts for 
male killing (MK). With probability (I — /3)(1 — r) create a new pair (to, /). 

(3) TOu, X /u, : same as above. 

(4) rriw X /: with probability q, create no offspring. This is the effect of cytoplasmic 
incompatibility (Cl). With probability 1 — q, create a new pair (to, /). 


Notice that, in contrast to m, the sex ratio at birth will not be I : 1, it is distorted 
due to male killing. Also we allow different mortality rates for males and females, in 
general. Therefore, even in the case when there is no male killing, the sex ratio would 
be distorted, in general. Also, we assume that any offspring resulting from Cl crossing 
is uninfected. 

We apply the mating rules described above to construct the birth function in our 
model. If the population sizes in the four compartments are denoted by M, , F, F ^, 
respectively, then the total number of possible matings is (M -|- Myj){F + Fyj). The total 
number of matings for example between uninfected males and infected females is MFw 
Hence the probability that a given mating of type mx takes place is ^ . 

To compute the total number of matings per unit time we follow the harmonic mean birth 
function approach from [22) . Accordingly, the total number of matings is proportional 
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M 


F + F,i. 


F 


M + Mjji] + -F + Fq^ 

hi hdqqq 

hi + hdqqj + F + Fqj^ 


hl„ 


F + F„; 


F, 


hi + hlqq, + F + Fqq 
hi + hlqq, 

hi + hdqjj + F + Fqq_ 


Jhd + hl„){F + Fq,) 

hi + hlqqj + F + Fqqq 

( 2 . 1 ) 


Hence the birth rate of offspring arising for example from the mating between an unin¬ 
fected male and an infected female is proportional to ■ We also naturally 

assume that there is competition between female individuals for finding an appropriate 
water reservoir to lay eggs. This is taken into account via a function \{Ftotai) which we 
assume (at least in the first instance) to be a monotonically decreasing function of the 
total number of females Ftotai, to allow for this competition for nesting places. Though 
\{Ftotai) is decreasing, it may approach a positive limit as Ftotai oo. This is to allow 
for the fact that gravid females that cannot find a place to lay their eggs may destroy eggs 
previously laid by others, and lay theirs instead. Thus the overall egg-laying rate should 
approach a positive limit as Ftotai —> oo, and therefore we assume that X{Ftotai) is a 
decreasing function such that A(c») > 0. Based on the individual mating rules explained 
earlier, our model reads as follows. 


hl\t) 

F'(t) 

hiut) 

Kit) 


-MmM + + (1 - /3)(1 - r)(MF™ + M^F^) + (1 - g)M„F), 


-HfF + 

f^mwhlq 


N 

XiFtotal ) 

N 

XtyFtotal) 


( 2 . 2 ) 


N 

X^Ftotal ) 


(MF + (1 - /?)(! - r)(MF„ + M„F^) + (1 - g)M„F) 
(1-/3 )t(1-7)(MF^ + M^F„), 


— fJ-fwFqjU H-- (1 — j3)T{hIFqqq + Mqq, Fqjj) . 


(2.3) 

(2.4) 

(2.5) 


A complete list of the variables, parameters and coefficient functions appearing in model 
(2.2)-(2.5) is given below. 


• M : number of uninfected male mosquitoes. 

• F : number of uninfected female mosquitoes. 

• h/Iqq, : number of Wolbachia infected male mosquitoes. 

• Fqqj'. number of Wolbachia infected female mosquitoes. 

• hi total = M + Alqjj, total number of male mosquitoes. 

• Ftotai = F + Fqqq, total number of female mosquitoes. 

• N = hitotai + Ftotai, total number of mosquitoes. 

• /3: reduction in reproductive output of Wolbachia infected females. 

• r: maternal transmission probability for Wolbachia infection. 
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• q: probability of cytoplasmic incompatibility (Cl). 

• 7 : probability of male killing (MK) induced by Wolbachia infection. 

• X{Ftotai)- average egg laying rate, which depends on the total number of female 
mosquitoes. 

• per-capita mortality rate of uninfected male mosquitoes. 

• fif: per-capita mortality rate of uninfected female mosquitoes. 

• per-capita mortality rate of Wolbachia infected male mosquitoes. 

• fifw- per-capita mortality rate of Wolbachia infected female mosquitoes. 


Model (2.2)-(2.5) is our starting point for a study of the Wolbachia infection dynamics in 
a sex-structured mosquito population. Later, we will expand this model by introducing 
WNv infection. 


2.2. Positivity and boundedness. First we begin by establishing positivity and bound¬ 
edness of solutions of model (2.2)-(2.5). 

Proposition 2.1. Assume that X is a monotone decreasing function such that 

lim Xi^FiQidi) — XjYiin ^ 0 , '^( 0 ) ^ 5 }7 ^ (^*^) 

Etotai-^OO 


hold. Then, the variables (M,F,Fu,) satisfying equations (2.2)-(2.5) remain non¬ 
negative if they are non-negative initially, and they remain bounded for all times. 

Proof. First note that the solution variables remain non-negative for all time; this follows 
from results in [53]. Adding equations (2.3) and (2.5), and noticing that /3 € [Oil]; 
q € [ 0 , 1 ], we have 

Hotai < - min{^/,/r/i„}Ftotai + - MtotaiFtotai 


Mi 


total 


total 


Therefore, 


< - mm{fJ.f,fif^}Ftotal + MtotaiFtotai 

^total 

— ( minl^xy, g.fw } F Xi^FiQi^i )) Fiqiq^i. 


lim sup Ftotai {t) < F 


(2.7) 


where F is such that A(F) = min{/r/, Note F exists since we assumed A is 

monotone decreasing and satisfies ( 2 . 6 ). 

Since Ftotai remains bounded it follows that Mtotai is bounded as well, because adding 


(2.2) and (2.4) we have 


Mtotai < - ^VCL{h-m, h^mw}Mtotai + MtotaiFtotai 

^yhotal + J^total 

^ “h F^Qidi 


( 2 . 8 ) 
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where B is an upper bound for X{Ftotai{t))Ftotai{t), which exists since Ftotai{t) is bounded 
and therefore so is X{Ftotai{t))Ftotai{t). From the differential inequality (2.8|, we can 
conclude that Mtotai is bounded, too. □ 

2.3. Boundary equilibria and their stability. It is straightforward to see that model 
(2.2)-(2.5) has only one non-trivial IFoZfeac/iza free boundary equilibrium if* = (M*,F*,0,0), 
where F* satisfies 

X{F*) = fif + n^, (2.9) 

and 


M* = 


hf* 

l^m 


( 2 . 10 ) 


under the assumptions that A(0) > /i/ + jim, A is a decreasing non-negative function, 
and X{Ftotai) -t Xmm (with X^in sufficiently small) as Ftotai -t oo. 

Note that there is no Wolbachia infected boundary equilibrium unless r = 1, a case 
that we shall treat separately later. 

Theorem 2.1. Suppose that X is a monotone decreasing non-negative function such that 
X{Ftotai) -t Xmin as Ftotai OO, with Xmin Sufficiently small, A(0) > pif + pim and 

1^/(1-^)'r 


< 1 . 


f^fw 


( 2 . 11 ) 


Then, the Wolbachia free boundary equilibrium E* = {M*, F*,0,0) of model (2.2)-(2.5) 
is locally asymptotically stable. 


Proof. Linearisation of system (2.2)-(2.5) at the equilibrium E* yields the following par¬ 
tially decoupled systems. The first system below is just the linearisation of equations 
(2.4)-(2.5) at the steady state, which we shall use to show that {My,{t), Eyj{t)) -A (0,0) 


as t 


oo. System ( |2.13[ ) is just system (2.2)-(2.3) in the case = Fy, = 0. 

-PmwMy, + - ^)^(1 - lW*Fy„ 



M* -\-F* 

—pfwFw + (1 ~ I3)tM Fy,, 

M' = -yimM 


F' = -pfF 

From the second equation of ( 2.12[ ), it is clear that if 

X{F*) 


M-hF 

^MF. 
M + F 


M* -h F 
oo. Then My,(t) 




( 2 . 12 ) 


(2.13) 


(2.14) 


0 as t —>■ 00 follows from the first equation of 


then Fy, (t) —)• 0 as t 

(2.12). Since M* and F* are given by (2.9) and (2.10), inequality (2.14) is equivalent to 
assumption (2.111. 
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It remains to prove the local stability of {M,F) = {M*,F*) as a solution of system 


(2.13). The Jacobian matrix of system (2.13) evaluated at {M*,F*) is given by 
j{M*,F*)= —^ 

/i/ + \ Mm fJ-ftJ-m + ^ {F )F fj,f 

The eigenvalues A of J {M*, F* ) satisfy the characteristic equation 

- X'{F*)F*fif)A - y{F*)F*{^lf + /rm)M/Mm = 0. 

Since A(-) is a non-negative decreasing function, X'{F*) < 0. We have 

Ai + A2 = - X'iF*)F*^if) < 0, 

and 

A1A2 = -X'{F*)F*{fif + Aim)M/Mm > 0, 

which implies ReAi < 0 and ReA2 < 0, so that {M*,F*) is locally stable as a solu¬ 


tion of (2.13). Therefore, the Wolbachia free equilibrium E* = (M*,F*,0,0) is locally 
asymptotically stable as a solution of the full system (2.2)-(|2.5|). □ 


If T = 1, i.e. we have complete maternal transmission of Wolbachia, then a boundary 
equilibrium of the form (0,0, M^,F^) may exist. The components of such an equilibrium 
solution must satisfy 


fJ-mwMw — (1 'l)h‘fwF,^, 


l-^fw — 


KF*) 


M* 


■F: 




(2.15) 


Moreover, (1 — 7)^/™ + fJ^mw = (1 — /3)(1 — a)HFw) must hold. Next we study the linear 
stability of such equilibrium, showing that it is linearly stable under condition (2.16) 
below. Inequality ( 2.16[ ) does not depend on 7, the male killing rate, but the steady state 
components and F^ do depend on 7 in the manner expected (for example, = 0 
when 7 = 1). Although Theorem |2.2| only apples if r = 1, we will be interested later on in 
the case when r is just slightly less than 1. Then, maternal transmission is imperfect and 
Wolbachia infected females produce small numbers of uninfected offspring. We anticipate 
that as T decreases from 1 to a value just less than 1, the equilibrium (0,0, M^,F^) shifts 
to another nearby position with small numbers of Wolbachia uninfected individuals and 
large numbers of infected ones; with no change of stability for r close enough to 1. 
The existence and stability of such an equilibrium will be important later on when we 
introduce West Nile virus (WNv) disease dynamics because, at a WNv-free equilibrium 
with large numbers of Wolbachia infected mosquitoes, the basic reproduction number 
Rq for WNv is likely to be less than 1. The implication is that Wolbachia infection in 
mosquitoes has the potential to control WNv infection. It does so by disrupting WNv 
virus replication causing WNv infected mosquitoes effectively to remain permanently (or 
for a very long time) in the latent stage of WNv. 
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Theorem 2.2. Assume that t = 1, \ is monotone decreasing with lim X{Ftotai) = 

Ftotal —>-00 

^min, (^min Sufficiently small) and 


A(0)> 


(1 - 'Y)Hfw + h-r. 


(1-/3)(1-7) 

holds. Then, an equilibrium of the form (M, F, M.^], Fw) = (0, 0, exists, and it 

is locally stable as a solution of (2.2) -([2^ if 

< {1- (2.16) 

Proof. The proof is similar to that of Theorem |2.1| The linearisation around (0,0, Mf ,, F ))) 
yields a system of linear equations for (M, F) that (when r = 1) are decoupled from the 
rest of the system. Moreover, it may be shown that (M, F) -A (0,0) as t —>■ oo if 

(1 - <7)M;A(F:) 


M* 


. p* 


< Tf 


holds, which becomes inequality (2.16) when the equilibrium equations (2.15) are invoked. 
Then, the and equations are considered, in the case when t = 1 and F = M = 0. 
Tedious computations yield that the linearisation of that system around the steady state 
has the Jacobian matrix equal to ((1 — 'y)fJ,fw + times 


(1 


f^mw 


(1 - /3)(1 - 7) [(l5^) l4^ + A'(F:)F:(1 - 

~ti-fwtimw + (1 ~ /3)A (F^,)F(„(1 — 7 )m/ui J 

and it may be further shown that its eigenvalues both have negative real parts. Thus we 
conclude that the steady state {0,0, Mf,, Ff,) is locally asymptotically stable. □ 

Next we prove that, under certain conditions, both infected and uninfected mosquitoes 
die out. Note, however, that (0,0,0,0) is not technically an equilibrium of (2.2)-(2.5). 

Theorem 2.3. Suppose that A is monotone decreasing and that 0 < X{Ftotai) < P-f + Tm 
for all Ftotai > 0, and that 

(Pf + Pni){^ - P)t < Pfw (2.17) 

Then {M{t),F{t),M.u,{t),Fy]{t)) -A (0,0, 0,0) as t ^ oo if all of the four variables are 
sufficiently small initially. 


(2.18) 

(2.19) 


Proof. From (2.4) and (2.5), 

Miit) < -pmwM.ui{t) + {pf + Mm)(l - P)r{l - ^)Fffit), 

F)„{t) < -pfwFw{t) + {pf + Pm){l - ffirFffit). 

From (2.19) and (2.17) we have Fy,{t) —>■ 0 as t > oo. Then inequality (2.18) implies 
that Mffit) -A 0 also. With = 0, we are reduced to system (2.13) and we 

now show that {M{t),F{t)) -A (0,0) as t —)■ oo, though this result is local, i.e. for small 
introductions of F and M. Note that {M,F) = (0,0) is not an equilibrium of (2.13), due 
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to the singularity, but we can remove that singularity by introducing the new variable 
^ = F/M. In terms of the variables ^ and M, system ( |2.13 1 becomes 

'1-CW' 


M'{t) = 


rnmimt)), 


( 2 . 20 ) 


We now show that (^, M) = (^*,0) is a locally stable steady state of (2.20), where 

_ -^(0) "f Mm ~ ^^f 

A(0) + ^/ - Aim’ 

provided > 0. The latter is not automatic. However, note that, from the first equation 


( 2 . 21 ) 


of (2.13l, M'{t) < —^rn)M{t). Therefore, if A(0) < /im then M{t) —>• 0. The second 

of (2.13) then gives F'{t) < —^fF(t) + \{Q)M{t) so that F{t) —>■ 0. By similar reasoning 
we arrive at the same conclusion if A(0) < Therefore, we may assume henceforth 
that A(0) > max(/i/. Aim) and, under these circumstances, > 0. The linearisation of 
the second equation of (2.20) near the equilibrium (^, M) = (^*,0) reads 


M'{t) = + A(0) j^^M(t) = i(A(0) - Aim - 

and therefore, since A(0) < Aim + M/i have M{t) —)■ 0. In this limit the ^ equation 
becomes 

eit) = -(/!/ - /im)?(t) + mm = mt)), 

and to show that is locally stable as a solution of this equation, it suffices to show 
that F'{^*) < 0. But, after some algebra, we have 

(^(0))^ - (/im -/i/)^ 

= -2AW-■ 

To show that F'{^*) < 0 holds, it suffices to show that A(0) > |Aim — fJ-fl, be. that both 
A(0) > Aim ~ Ai/ and A(0) > /i/ ~ fJ-m hold. But this follows from the fact that we are now 
restricting to the case when A(0) > max(At/, Aim)- Therefore, the proof of the theorem is 
complete. □ 

2.4. Existence of strictly positive steady states. In this section we examine the 
possible existence of coexistence steady states {M,F,M^,F^) of model (2.2)-(2.5), i.e. 


steady states in which each component is strictly positive. It turns out that in some 
parameter regimes multiple coexistence steady states may exist while, in others, there 
is just one or none at all. An understanding of these properties helps us to understand 
how one might exploit Wolbachia infection in mosquitoes to effectively control WNv. In 
this section we simplify by assuming that 7 = 0, i.e. that there is no male killing. 
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At the steady state, dividing ( |2.2[ ) by ( |2.3| ), and \2A\ by ( |2.5[ ), we obtain 

. ib 


M = F^, 


l^T) 


( 2 . 22 ) 


From (2.3) and (2.5) we then obtain 
M/F = 


X + {^FF^ + ^Fl\+{l-q) 

\ l^m \ l^m l^mw / 


fj-fu 


-ff, 






l-^m l^mw 


(2.23) 

(2.24) 


respectively, where A* = X{F + Fyj). From ( |2.24[ ) we have 


Fyj I fifw A 






+ F I IJ-fw + 


Mr) 


= A* ( F(1 - /3)t^ + F^(l - f3)T^ 




From (2.25) we obtain 
where 

ki(A*) = fj.fw + - A*(l - 

l^mw l^mw 


FyjKl{K) = Fk2{X*), 


(2.25) 

(2.26) 


LL f LL f 

/^2('^*) - l^fw l^fw • 

l^m 


(2.27) 

Note that if at the strictly positive steady state, we have = 0, then this necessarily 

implies that ^s:2(A:^) = 0. This is only possible if 

f^fw l^mw 

M/ Mm 

holds. This case is excluded from Theorem 12.41 but is treated in the next subsection. 


(2.28) 


If (2.28) does not hold then, using (2.27), from (2.23) we obtain 

l^fw 




0 =Ki(A*) I fj,f + — - A*— j - K2(A*)A*(1 - /3)(1 - r)- 

\ Mm Mm I t^mw 

+ ki(A*)k2(A*) (^J.f + - A*(l - ^)(1 - t)^ - A*(l - 


(2.29) 

l^mw f-^m l^mw / 

The right hand side of ( |2.29 1 is, in general, a cubic polynomial in A*. If there exists a 
positive root Aj, then since A is a strictly monotone function, a corresponding unique 
F^ + F^ value may be found. From (2.23) we may then determine a unique solution 
{F"^,F^). Further analytic progress is possible in certain particular cases of interest, 
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which we now investigate. The first concerns the case when r = 1, or when r is very close 
to 1, meaning that maternal transmission of Wolbachia is complete or nearly complete. 
This is in fact the biologically relevant case for a number of Cl inducing Wolbachia 
strains in mosquito species treated in the literature, see e.g. [9]. This leads us to expect 
the existence of a steady state of (2.2)-(2.5) with large numbers of Wolbachia infected 
mosquitoes and few, or no, uninfected ones. The stability of such a steady state still 
depends on the other parameter values, and it will be stable if there is a high probability 
of mating between infected males and uninfected females resulting in no offspring (the 
effect of Cl), i.e. q is close to 1; see also inequality (2.16) for the case r = 1. The 


existence of a stable steady state of (2.2)-(2.5) with the above mentioned properties 


is important because the low number of Wolbachia uninfected females implies that the 
quantity F*, featuring in the first term of the parameter i?o defined later in (3.47), 
is small. The likelihood of i?o being less than 1 (the condition for WNv-eradication) 
depends mostly on that first term involving F* ^ since the second term in (3.47) involves 


a small parameter e and is automatically small. For these reasons, we are interested in 
stable steady states of model (2.2)-(2.5) of the form (M*, F 


M^,F^), with M* and F* 


small compared to and F^ (and, ideally, M* = F* = 0). Therefore, referring to 
Theorems |2.1| and |2.2[ and restricting for now to the case r = 1, we ideally would like the 
boundary equilibrium (M*,F*,0,0) of Theorem 2.1 to be unstable, and the boundary 


equilibrium {0,0, M^, F^) of Theorem 2.2 to be linearly stable. The conditions for this, 
when T = 1, are that the inequalities ^/(l —/3) > /i/u, and {l — q)fifw < {1 — /3)fif should 
hold simultaneously, but note that the second of these follows from the first. For r = 1, 
guided by elementary competition theory, instability of one boundary equilibrium and 
stability of the other suggests that there will be no coexistence equilibrium, and this is 
what we prove in Theorem |2.4| below. If r is decreased from 1 to a value slightly below 
1, there is no longer a boundary equilibrium with only Wolbachia infected mosquitoes 
present. What happens is that the equilibrium {0,0, Mi^, F^), which exists when t = 1, 
moves to another nearby point in so that Wolbachia infected mosquitoes now coexist 
with uninfected ones, the former being dominant. Very importantly, this will be the only 
coexistence steady state if r is sufficiently close to 1, and it is a desirable steady state 
for WNv eradication because particular Wolbachia strains can significantly reduce WNv 
virus replication in mosquitoes, see [18]. 

Theorem 2.4. Suppose that 7 = 0, t = 1, /r/(l — /3) > p-fw, PmPfw 7^ PfPmw, and 
that X is a strictly positive decreasing function with A(oo) = Xmin (and Xmin sufficiently 
small). Then system (2.2)-(2.5) has no coexistence equilibrium with M*, F*, M)), Ff > 0. 


If the foregoing hypotheses hold, except that t is slightly less than 1, then system (2.2)- 


(2.5) has precisely one coexistence equilibrium in which Wolbachia uninfected mosquitoes 


exist in very small numbers relative to Wolbachia infected ones. 


Proof. Since we assume r = 1, the form of (2.29) simplifies and in fact we may cancel 


ki(A*) since we seek equilibria in which M*, F*, M*, F* > 0. After some further algebra. 
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we find that there is just one value for A*, given by 


l^mw \ l^mw /^m / ^ml^mw 

I M/M/ii) _ 




(2.30) 


l-^m l^mw 

If A* < 0 then the equation A* = X{F + F^,) cannot be solved for F + Fu,, so we may 
restrict to the case that A* > 0. Recalling that ki and K 2 are defined by (2.27), we find, 
with A* given by (2.30) and with r = 1, that 

\ 1 


Ki(A*) = 1 - 




^■fw + ^^fw) 


and 


K2(A.) = i((l-«w-f.,0 

Q \ / 


(2.31) 


(2.32) 


Since /r/(l — l3) > /i/u, it follows that the sign of the product ki(A*)k 2 (A*) is the same 
as the sign of 


1 - 




1 - 






and, since we assume ^ P^fl^mw, it follows that ki(A*)«; 2 (A*) < 0. This makes 

it impossible to find F > 0 and F^, > 0 satisfying (2.26), and so there is no coexistence 
equilibrium. 

Next we prove the second assertion of the theorem. Let r = 1 — e. Since we expect the 
equilibrium (0, 0, M)),, F)^), which exists when r = 1, to move to another nearby point 
when r = 1 — e, for e sufficiently small, we seek an equilibrium of (2.2)-(2.5) of the form 

M(e) + • • • , 

F(e) =eF^^'> + + • • • , 

M^{e) =M; + + ... , 

Fwis) =F^ + sFj;^^ + e'^F^^ H-. 

Coefficients of £ yield 


and 


Using 


= 


HF^) 

MZ + F* 1 

MFZ) 

Ml + F* 


(1 - /3)m;f: + (1 - (7)m;f(i) 
(1 - I3)M:f: + (1 - g)M:F(i) 


(2.33) 

(2.34) 


„ A(F*) 

~ M* ' ^ P)^w^ 
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(2.34) reads 


so that 


fifl 




^( 1 ) = 


1 -, 

(1 - 


(2.35) 


Note that F^^"> > 0 because of the assumption /i/(l — P) > hence from (2.33) we 
conclude that > 0 holds, and 

— P)F^ 


= 


(2.36) 


Mm [(1 - /3 )m/ - (1 - q)^J-}w\' 

In conclusion, if r is reduced from 1 to the value 1 — e then the equilibrium (0,0, M*, F*) 
moves to {eM^^\eF^^\ + eMw \ F^ + eF^'^), with M^^'> and F^^'> given by (2.36) 

and (2.35[), respectively. □ 


Note that from the proof of Theorem 2.4 we can see that at the equilibrium the 


male/female ratio for Wolbachia uninfected mosquitoes is given approximately by 


Rl_ 


which we treat now. In 


In this situation, we have 


Theorem |2.4| excludes the case when /imM/u; = M/Mm 
particular, we assume that 7 = 0 , and — = fi = thisL 
M = fiF, and From equations ( |2.3[ ) and ( |2.5[ ) we obtain 

, zr 1 + (1 - /3)(1 - r) (FF^ + F^) + (1 - q)FF^) , (2.37) 


f^fw^w - 


(1 + /i)(F + Fuj) 


((l-/3)r(FF^ + F2)). 

(l + M)M/ii 


(2.38) 


(1 + /r)(F + Fyj) 

From equation ( 2.38| ) we find that 

A* = A(F + F^) = ^i±iffe, (2.39) 

M(I - P)T 

hence for the existence of a coexistence steady state it is necessary that A(0) > 77 ^ 73 ^^, 
in which case there exists a unique c = F + F^, such that A(c) = holds. Then, 

from equation ( 2.37|) we obtain 


M/ 

f^fw 


c(I - P)tF = F" + (1 - q)FF^ + c(I - /?)(! - r)F^, 


(2.40) 


from which, using F^ = c — F, we obtain 


0 = F^g + Fc f(l - g) - (I - /3)(1 - r) - (1 - P)t-^) + c^(l - /3)(1 - r). (2.41) 

V M/»/ 

For the existence of the positive steady state we need to guarantee that the quadratic 
equation above has (at least one) positive (real) solution, and that the solution is less 
than c. 
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From equation (2.411 it is clear that in case of complete maternal transmission, i.e. 


for T = 1 there cannot be more than one coexistence steady state. In this case, from 


equation (2.41) we obtain 


F = 




(2.42) 


Therefore, if 


l-q< (1-/3)^ < 1 


fj-fv 

holds, then we have 0 < F < c, and if A(0) > 
positive steady state exists. 


(2.43) 

also holds, then a unique strictly 


Circumstances under which (2.43) is likely to hold include that q is sufficiently close 


to 1 and, at the same time, /? is sufficiently close to 1 or is less than 1. The biological 
interpretation is clear: for the existence of a coexistence steady state, the fertility cost 
(or mortality increase) due to Wolbachia infection should be sufficiently large. 

It is clear from (2.41) that we can never have more than two coexistence steady states. 


On the other hand it is interesting to show that for some realistic parameter values it is 
possible to have two coexistence steady states. 

To this end we consider the case q — 1, = 1, and we assume that r ^ 1, ,5 ^ 1. 


In this case, from (2.41), we have 


Fi/2 = C- 


That is, for 0 < F 1/2 < c to hold, we need to assume that 




(2.44) 


(2.45) 


hold simultaneously. It is easy to verify that this can be achieved for any c, for example 
with r = 0.99, (3 = 0.5. Note that the condition = 1 can be relaxed, too. Also, by 
continuity arguments, it follows that for parameter values close enough we still have two 
coexistence steady states. We summarise our findings in the following proposition. 


Proposition 2.2. In the case when 7 = 0, g = 1 and — = there exists a set of 
values for the remaining parameters such that system (2.2)-(2.5) admits two coexistence 
steady states. 


In summary, we have shown that our Wolbachia model (2.2)-(2.5) may exhibit all of 


the three qualitatively different possible scenarios, i.e. when there are 0,1 or 2 coexistence 
steady states. 
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3. Model incorporating West Nile virus (WNv) 

There has been considerable recent interest in West Nile virus (WNv), with the great 
majority of mathematical papers on the topic having appeared in the last 15 years. Nu¬ 
merous types of models have appeared, some including spatial effects and others giving 
consideration to issues such as age-structure in hosts, optimal control or backward bifur¬ 
cation. See, for example, Blayneh et al [3], Bowman et al [1], Gourley et al El, Lewis 
et al [23j and Wonham and Lewis m- 


We introduce the following model as an extension of model (2.2)-(2.51 to include WNv 


dynamics. Our model for WNv dynamics has similarities to that in Bergsman et al [1] 
with one resident bird population. Hence, as in [1] and in some references therein, we 
compartmentalise the vector and bird population into SEI and SEIR classes, respectively. 
Taking into account our earlier model the complete WNv- Wolbachia mosquito-bird pop¬ 
ulation model takes the following form. 


F' = 




M Mw F Fi 

- ^ifFs - afPbfFs 

B 


-{MF + (1 - /3)(1 - t){MF^ + M^F^) + (1 - q)M^F) 

I 

B, 


K = o^fPbfFs- 
F'i = 


M + Muj -\- F -\- 


Btotal 
l^fFe ^fFe^ 

-(1 F^n + MyjF^^ l^fwFujS ^fwPbfF-u) 




Btotal 


Fyje ~ ^fwPbfF^ 


’ D 

^total 

Fyji f Fyje. PfwF'ujii 


PfwF'nje fFyjQ^ 


Bg — Yl{^Btotal') PbBs (CUfPfhFi “h OifwPfjjF^ 

\ Btotal 


B.. 


Btotal J 


H' = afPfbF^ 


B., 


^ f wP fbFwi 


B, 


I —J UJjr J U- UJL 

^total ^total 

B^ — l^bBe PbBi 
Bj. = i^iBi 

^{Ftotal^ 


- flbBe - nBe 


M' = 
M' = 


Af + Mu} F F F F-ii 

^{Ftotal') 

M F Myj F F F F-ii 


(MF + (1 - /5)(1 - r){MF^ F M^F^) + (1 - q)M^F) - 
(1 - /5)t(1 - j){MF^ F M^F^) - (3.46) 
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where 


Btotal — Bg+Be + Bi+Br, 

Ftotai = F -\- Fiu, Fjj, = F^s + F^f, + F^i, F = Fg + F^ + Fi. 


Some of the parameters of system (3.46) are defined after (2.2)-(2.5); the rest are defined 
as follows: 


• af'. biting rate of female Wolbachia uninfected mosquitoes; 

• a/u,: biting rate of female Wolbachia infected mosquitoes; 

• pbf'. transmission probability of WNv from infectious birds to WNv-susceptible 
female mosquitoes; 

• pfb- transmission probability of WNv from WNv-infectious female mosquitoes 
to susceptible birds; 

• lyf. per-capita transition rate of WNv-exposed female Wolbachia uninfected 
mosquitoes to the infectious stage of WNv; 

• e € [0,1]: small parameter modelling increased time that Wolbachia infected 
mosquitoes spend in the latent stage of WNv, due to the tendency of Wolbachia 
infection to hamper the replication of WNv in mosquitoes; 

• Vb- per-capita transition rate of WNv-exposed birds to the infectious stage of 
WNv; 

• lyp. per-capita rate at which infectious birds recover; 

• pb- per-capita natural death rate for birds; 

• Phi' per-capita WNv-induced death rate for infectious birds. 


In this section, it must be emphasized that we are considering two different kinds of 
infection. Mosquitoes may be infected by either Wolbachia or WNv, or both. WNv 
infection is assumed possible only for female mosquitoes (since it is females that bite) 
and is modelled using an SEI (susceptible-exposed-infectious) approach with subscripts 
s, e and i. The variables Fg, Fg and Fi denote the numbers of Wolbachia uninfected 
mosquitoes that have susceptible, exposed and infectious status with respect to WNv. A 
subscript w indicates Wolbachia infection, so that Fws^ F^g and F^i denote the numbers 
of Wolbachia infected mosquitoes that have susceptible, exposed and infectious status 
with respect to WNv. The variables M and are the numbers of Wolbachia uninfected 
and Wolbachia infected male mosquitoes, none of which have WNv. Birds are only 
susceptible to WNv and their numbers are given by the variables Bg, Bg, Bi and Br 
denoting susceptible, exposed, infectious and recovered birds. Many of the terms of 
system (3.46) are also present in (2.2)-(2.5) without change, here we just discuss the 
extra terms that model the addition of WNv dynamics. WNv-susceptible mosquitoes, 
whether Wolbachia infected or not, acquire WNv infection by biting infectious birds; this 
is modelled via the last term in the first and fourth equations of (3.46) using the idea of 
mass action normalised by total host density, the biting rates (the a parameters defined 
above) having been separated out, rather than being absorbed into the transmission 
coefficients pbf and pfb as is often customary. Having acquired WNv infection from a 
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bird, a mosquito enters the latent phase of WNv and is classed as an exposed mosquito. 
Exposed mosquitoes become WNv-infectious at rates I'fFf, and svfFy^e for Wolbachia 
uninfected and Wolbachia infected mosquitoes, respectively. In the latter, the presence 
of £ € [0,1] models the tendency of Wolbachia infected mosquitoes that have contracted 
WNv infection to spend a greater amount of time in the latent stage of WNv, since 
Wolbachia infection tends to block WNv replication making it less likely that such a 
mosquito would ever become WNv-infectious. Of course, we are at liberty to take e very 
small indeed, with the implication that the Wolbachia infected mosquito spends so long 
in the latent stage of WNv that it probably dies in that stage. This is our approach to 
modelling the blocking of WNv replication by Wolbachia. 

Birds acquire WNv from bites by WNv-infectious mosquitoes, which may or may not 
have Wolbachia infection as well. Thus there are two infection rates for birds, these 
can be found in the right hand side of the seventh equation of (3.46), and also in the 


eighth equation since birds initially enter the exposed stage of WNv. This has a mean 
duration of l/vt for birds, after which they become WNv-infectious. Birds may recover 
from WNv, at a per-capita rate Note that, for birds, death due to WNv is modelled 
using a separate parameter /if,i to distinguish from natural death, accounted for by p,f,. 
The function Il{Btotai) is the birth rate function for birds. 

The approach we use here to model the latency stage of WNv (in either birds or 
mosquitoes) is not the only possible approach. Our approach permits individuals to 
spend different amounts of time in the latency stage, and we may only speak of the 
mean time spent in that stage. There are other approaches in which all individuals of 
a particular status (for example, all Wolbachia uninfected mosquitoes) spend the same 
amount of time in the latent stage of WNv. The time could be different for Wolbachia 
infected mosquitoes. These approaches result in models with time delays. 


3.1. Local stability of the WNv-free equilibria. Equilibria of system (3.46) may ex¬ 
ist in which WNv is absent. Such WNv-free equilibria include the equilibrium {M*,F*,0, 0) 
considered in Theorem |2.1[ in which both WNv and Wolbachia are absent, and equilibria 
in which WNv is absent but Wolbachia are present. We show that multiple WNv- 
free equilibria may coexist that have both Wolbachia uninfected and Wolbachia infected 
mosquitoes, we present a necessary and sufficient condition for any particular WNv-free 
equilibrium to be locally stable, and we show that the most likely scenario for eradica¬ 
tion of WNv is to have large numbers of Wolbachia infected mosquitoes, with solutions of 


system (3.46) evolving to a WNv-free equilibrium that has large numbers of Wolbachia 


infected mosquitoes and relatively few uninfected ones. Theorem 3.1 applies to any 
WNv-free equilibrium, of which there may be several. Of course, we may have i?o < 1 
at one WNv-free equilibrium and Rq > 1 dX another. It depends on the values of F*, 
F^g and B* for the particular WNv-free equilibrium under consideration. For clarity 
of exposition, we include as a hypothesis that the equilibrium be stable to the subset 
of perturbations in which WNv is absent (i.e. stable as a solution of the subsystem 


(2.2)-(2.5)), rather than including explicit conditions for stability of an equilibrium as a 
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solution of that subsystem. The latter stability problem is a tedious one in its own right 
and is under consideration elsewhere in this paper. Theorem |3.H in the form presented 
below, highlights clearly the particular role played by Rq. 

Theorem 3.1. Let F*, and B* be the equilibrium values for the female susceptible 
(to WNv) Wolbachia uninfected and Wolbachia infected mosquitoes and susceptible birds, 
in any WNv-free equilibrium. Let 


Rq — 


_ I'bi'fPfbPbf _ / ctjjFs /Bt) 

iHb + fJ-bi + + >^b) + pLfwiPf'u 




(3.47) 


Then, if Rq < 1, the WNv-free equilibrium under consideration is locally stable as a 


solution of the full system (3.46), if it is stable to perturbations in which the exposed and 


infectious variables remain zero. 


Proof. At any WNv-free equilibrium, the linearisation of system (3.46) decouples to some 


extent making it sufficient to show that, when Rq < 1, each component of the solution 
of the following system: 

ajpbfFf 


K 

F( 

F' 

we 

FL 

B'e 

B' 


B\. 

^fwPbfF^ 


Fi i.Pf F 






^^fF'we Pf wFyji, 

(AfPfbFi “t“ (AfwPfbFwi F b^b^B^ 

b^bBe (.hb F P’bi F b^i^Bi, 


(3.48) 


tends to zero. It is taken as a hypothesis that the susceptible variables then approach 


their respective steady state values. Note that system (3.48) has a structure that allows 
the application of Theorem 5.5.1 in Smith |29j . making it possible to restrict attention to 


the real roots of the characteristic equation associated with (3.48). That characteristic 


equation, corresponding to trial solutions with temporal dependence exp(At), is most 
easily analysed when written in the form 

{X F fJ-b F ^bi F b>i)[X fXb F Vb)B* 


= ^bVfPfbPbf 


Pb 


prv'^ 

fw^ ws 


(A F /^/)(A F F b>f) (A F ^fw){X F fifw F evf) 


(3.49) 


As functions of the real variable A, the right hand side of (3.49) is decreasing, at least 


for A > 0, while the left hand side is a quadratic with two real negative roots. A simple 
graphical argument shows that if the left hand side exceeds the right hand side when 


A = 0 (i.e., if Rq < 1, with Rq defined by (3.47)), then any real roots of the characteristic 


equation are negative which, since only the real roots need to be considered, implies that 
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Symbol 

Definition 

Value 


Per-capita mortality rate of male mosquitoes 

1/20 

M/ 

Per-capita mortality rate of female mosquitoes 

1/20 


Per-capita mortality rate of FP-infected male mosquitoes 

1/20 

f^wf 

Per capita mortality rate of FP-infected female mosquitoes 

1/20 

r 

Maximum per-capita mosquito egg-laying rate 

30 

k 

Competition coefficient for mosquitoes 

5000 

P 

Fitness cost of FP-infection on reproduction 

G [0,1] 

T 

Maternal transmission rate of Wolbachia 

G [0,1] 

q 

Strength of Cl due to FP-infection 

G [0,1] 

7 

Male killing rate due to FP-infection 

G [0,1] 

Uf 

Per-capita FP-free mosquito biting rate 

0.09 


Per-capita FP-infected mosquito biting rate 

0.09 

Pbf 

WNv transmission coefficient from birds to mosquitoes 

0.16 

Pfb 

WNv transmission coefficient from mosquitoes to birds 

0.88 

Pb 

Natural per-capita mortality rate of birds 

1/(365 X 3) 

f^bi 

WNv-induced per-capita death rate of birds 

0.1 

n{B) 

Birth rate of birds 

100/365 


Per-capita rate at which FP-free mosquitoes 



complete WNv-latency and become WNv-infectious 

G [0,1] 

£ 

eiyf is the per-capita rate at which FP-infected mosquitoes 
complete WNv-latency and become WNv-infectious 


Vb 

Per-capita rate at which exposed birds become infectious 

0.2 

Vr 

Per capita rate at which infectious birds recover 

0.2 


Table 1. Definition of parameters. Time is measured in days, and rates 
in day“^. W stands for Wolbachia so that, for example, IT-infection 
means Wolbachia infection. Parameter values have been chosen purely 
to demonstrate possible solution behaviour and are not based on data. 


each component of the solution of (3.48) approaches zero as t 
theorem is now complete. 


oo. The proof of the 
□ 


4. Numerical simulations 

In the simulations we set \{Ftotai) '■= where r is the maximum per- 

capita mosquito egg-laying rate, and k measures intra-specific competition among female 
mosquitoes. It should be noted that our model assumes that the mosquito population 
persists annually, as for example in the tropical climates of South-East Asia. 






N of Mosquitoes at steady state N of Mosquitoes at steady state 
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Figure 1. The dependence of the values of F and at the steady 
state of model (2.2)-(2.5) on the parameters /3, t , 7 and q. Here, all the 
per-capita mortality rates of mosquitoes were taken as 0.05, and we set 
r = 30, /c = 5000. 
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Figure 2. Simulation of model (2.2)-(2.5l showing F and against 
time. Here, all the per-capita mortality rates of mosquitoes are taken as 
0.05, and we set r = 30, fc = 5000, /3 = O.l, q — 0.9, 7 = 0 and r = 0.95. 



£ 


Figure 3. The basic reproduction number Rq, defined in (3.471, plotted 
against e. The parameter values are (3 = 0.1, q = 0.9, 7 = 0, r = 0.9, 
with the other parameter values given in Table For these parameter 
values, F* = 3630, = 23837, B* = 300, and almost all of the 

mosquitoes are infected with Wolbachia. In this case, if £ < 0.03, the 
basic reproduction number i?o < 1, and the WNv will die out. 















N of WNv-infectious, W-uninfected female M F? ^ susceptible birds B 
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(a) Bs versus time (b) Bi versus time 




(c) Fi versus time 


(d) Fyji versus time 


Figure 4. Simulation of model ( |3.46 ) 
given in the caption of Figure for the 
the case e = 0.02, WNv dies out. 


with the parameter values as 
cases e = 0.2 and e = 0.02. In 
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T 


Figure 5. Basic reproduction number Rq, plotted against r, with P = 
0.1, q = 0.9, 7 = 0, e = 0.02, and the other parameter values given in 
Table [2 If t > 0.9 the basic reproduction number Rq < 1, and WNv 
will die out. 



q 


Figure 6. Basic reproduction number Rq, plotted against q, with /3 = 
0.1, T = 0.9, 7 = 0, e = 0.02, and the other parameter values given in 
Table If g > 0.8 the basic reproduction number Rq < 1, and WNv 
will die out. 
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Figure 7. Basic reproduction number i?o, plotted against P, with r = 
0.9, q = 0.9, 7 = 0, e = 0.02, and the other parameter values given in 
Table [2 If /3 < 0.3 the basic reproduction number Rq < 1, and WNv 
will die out. 



T 

Figure 8 . Basic reproduction number plotted against 7 , with /3 = 
0.1, q = 0.9, r = 0.9, e = 0.02, and the other parameter values given in 
Table [2 If 7 < 0.65 the basic reproduction number Rq < 1, and WNv 
will die out. 


5. Conclusion 

In this paper we have derived a detailed sex-structured model for a mosquito popula¬ 
tion infected with Wolbachia . The model captures many of the well-known key effects of 
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Wolbachia infection, including cytoplasmic incompatibility, male killing, reduction in re¬ 
productive output and incomplete maternal transmission of the Wolbachia infection. Our 
analysis shows that the mosquito population can stabilise at a Wolbachia free equilibrium 
under certain circumstances, which include situations when inequality (2.111 holds. Such 
circumstances include, for example, if Wolbachia infection significantly reduces reproduc¬ 
tive output, and/or Wolbachia infection significantly lowers female life expectancy. We 
also showed that if t = 1, i.e. maternal transmission of Wolbachia is complete, then the 
mosquito population can stabilise at an equilibrium in which all mosquitoes are infected 
with Wolbachia. This happens in the case of sufficiently high cytoplasmic incompatibil¬ 
ity. In the case of r close to 1 we have shown that Wolbachia infected mosquitoes can 
coexist with small numbers of uninfected mosquitoes. We have also shown that under 
some additional assumptions our model has multiple coexistence steady states. 

We extended the sex-structured mosquito population model (2.2)-(2.51 to include West 
Nile virus, which is spread by birds and mosquitoes, treating WNv as an SEI infection for 
mosquitoes, and as an SEIR infection for birds. We were motivated by results recently 
reported in [TS] , which suggest that a particular strain of Wolbachia substantially reduces 
WNv replication in the mosquito species Aedes aegypti. We modelled this crucial phe¬ 
nomenon by incorporating a small parameter e, the reciprocal of which is proportional 
to the time spent in the WNv exposed class for Wolbachia infected mosquitoes. This 
enabled us to assess the potential of Wolbachia infection to eradicate WNv via its effect 
on WNv replication in Wolbachia infected mosquitoes. Notably the expression we ob¬ 
tained for the basic reproduction number Rq suggests that WNv will be eradicated if at 
the steady state the overwhelming majority of mosquitoes are infected with Wolbachia., 
and the Wolbachia infection substantially reduces WNv replication in mosquitoes. The 
first of these hypotheses is in fact shown to hold for a number of Wolbachia strains and 
mosquito species, see e.g. [8]. 
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